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Abstract. We obtain abundances of a, iron peak and neutron capture (n-capture) process elements in four Ba 
Q ■ stars HD 26886, HD 27271, HD 50082 and HD 98839 based on high resolution, high signal-to-noise spectra. We 

find that all of these Ba stars are disk stars. Their a and iron peak elements are similar to the solar abundances. 
The n-capture process elements are overabundant relative to the Sun. In particular, the second peak slow neutron 
capture process (s-process) elements, Ba and La, are higher than the first peak s-process elements, Y and Zr. 
Analyzing the abundances of four sample stars, the heavy-element abundances of the strong Ba star HD 50082 
^ ' are higher than those of other three mild Ba stars. The stellar mass of the strong Ba star HD 50082 is 1.32M0 

^S) I (-1-0.28,— O.22M0), which is consistent with the average mass of strong Ba stars (I.5M0). For mild Ba star 

. HD 27271, we derive I.9OM0 (-|-O.25,-O.2OM0), consistent with the average mass of mild Ba stars (1.9M0, with 

' O.6M0 white dwarf companion). For mild Ba star HD 26886, the derived 2.78M0 (-|-O.75,-O.78M0) is consistent 

with the average 2.3M0 of mild Ba stars with O.67M0 companion white dwarfs within the errors. Mass of mild 
Ba star HD 98839 is high to 3.62M0, which inspires more thoughts on the formation of Ba star phenomenon. 
Using our angular momentum conservation theoretical model of wind accretion of Ba binary systems, we obtain 
the theoretical heavy-element abundances of Ba stars that best fit our data. The results show that the observed 
• abundances of the typical strong Ba star HD 50082 and the typical mild Ba star HD 27271 are consistent with 

the theoretical results very well. This suggests that their heavy-element abundances were caused by accreting the 
Q . ejecta of AGB stars, the progenitors of the present white dwarf companions, through stellar wind. However, wind 

' accretion scenario cannot explain the observed abundance pattern of the mild Ba star HD 26886 with shorter 

I orbital period (P=1263.2 d). The mild Ba star HD 98839 with high mass (up to 3.62M0) and very long orbital 

. period {P >11000 d) may be either a star with the heavy elements enriched by itself or a "true Ba" star. 
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1. Introduction and subsequently mixed to the surface (e.g. Burbidge & 

Burbidge 1957 for HD 46407). However, the absence of the 

The "barium stars", or "Ban stars" as they were first , , ggrri / o., ms ^ 1 -.j 

„ , , „ , , „ ,. . , unstable nucleus ^''ic (ti = 2 x 10° yrs) supply evidence 

' ^ & y that they are not presently undergoing nucleosynthesis ac- 

Bidelman & Keenan (1951). These G and K giants ex- , , , .1-11 1 ri j- 

^ ' ° tivity and the third dredge-up process like asymptotic gi- 

hibit enhanced features of Ban, Sril, CH, CN, and some- branch (AGB) stars 
times C2 lines. Many qualitative studies on Ban stars 

have been developed (e.g. Garstang 1952; Burbidge & ^^er Burbidge et al. (1957) suggested that elements 

Burbidge 1957; Danziger 1965; Pilachowski 1977; Tomkin 'heavier than iron are synthesized in the interior of AGB 

& Lambert 1983; Smith 1984; Kovacs 1985; Zacs 1994). stars, it is generally beheved that Ba stars belong to bi- 

The results confirmed that, in addition to Ba and Sr, other ^ary systems and their heavy-element overabundances are 

heavy elements are also enhanced, e.g. Y, Zr, La, Ce, Pr, Produced by accreting the matter ejected by the com- 

Nd and Sm panions (the former AGB stars, now evolved into white 

Some researchers suggested that the overabundances ^^^^^f^)' ^^^^^ researchers have studied the binarity or 

of Ba stars were produced in the interior of the Ban star, heavy-element abundances of Ba stars (McClure et al. 

1980; McClure 1983; McClure & Woodsworth 1990; Boffin 

Send offprint requests to: Y. C. Liang, email: & Jorissen 1988; Jorissen & Mayor 1992; Jorissen et al. 

Yanchun.Liang@obspm.fr, lyc@yac.bao.ac.cn 1998; Liang et al. 2000). The mass exchange took place 
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about 1 X 10^ years ago, so the ^^Tc produced in the orig- 
inal thermal pulse (TP) AGB stars have decayed. The 
accretion may be in the forms of wind accretion, disk 
accretion or common envelope ejection (Han et al. 1995; 
Jorissen et al. 1998; Liang et al. 2000). 

Therefore, it is very important to combine abundance 
analysis and binary characters of Ba stars to study their 
properties and formation scenarios. 

At present, there is a large sample of Ba stars with 
measurements of orbital elements (Carquillat et al. 1998; 
Udry et al. 1998a, 1998b; Jorissen et al. 1998), abso- 
lute magnitudes and kinematics (Gomez et al. 1997; 
Mennessier et al. 1997). However, the corresponding 
heavy-element abundances have not been obtained from 
observations. To explore their formation scenario, it is nec- 
essary to study their observed abundances, combined with 
their orbital elements. 

Moreover, using our angular momentum conservation 
model of wind accretion of Ba stars formation, and com- 
bining the AGB stars nucleosynthesis (Liang et al. 2000; 
Liu et al. 2000), we can calculate the theoretical heavy- 
element abundances of Ba stars. Thus we can understand 
the observed abundance patterns of the sample stars from 
theory. 

In addition, with the high precision Hipparcos data, 
the precise photometric parameters, improved methods 
to determine stellar atmospheric parameters and devel- 
oped stellar evolutionary tracks study, it is possible to 
allow to obtain the masses and atmospheric parameters 
of these stars. For Ba stars, combining masses, element 
abundances and orbital parameters provides a good under- 
standing of their properties and formation. In this paper, 
we discuss the properties of four Ba stars by combining 
these parameters. 

This paper is organized as follows. In Sect. 2-4, we 
describe the observations and analysis methods in de- 
tails, and present the derived atmosphere parameters and 
masses. Detailed spectral lines analysis and equivalent 
widths (EWs) are also presented. In Sect. 5, wc show the 
element abundance patterns and analyze the abundance 
uncertainties. In Sect. 6, we calculate the theoretical abun- 
dance trends of s-process elements in Ba stars, and com- 
pare them with the observed patterns of our sample stars. 
Detailed discussions are given in Sect. 7. In Sect. 8, we 
summarize this paper and describe possible future studies. 

2. Observations and data reduction 

The high resolution spectra of four barium stars were 
obtained with the Coude Echelle Spectrograph and a 
1024x1024 TeK CCD attached to the 2.16m telescope 
at the National Astronomical Observatories (Xinglong, 
China) in 2001. The red arm of the spectrograph with 
a 31.6 grooves/mm grating was used. With a 0.5 mm slit, 
the resolving power is of the order of 40 000 in the middle 
focus camera system. The total wavelength coverage was 
5800-9000 A . The signal-to-noise (S/N) ratio is greater 
than 150 over the whole region. A detailed description of 




6140 6142 6144 

Wovelength (ongstrom) 



Fig. 1. Example of spectra of HD 98839: Ball A6141.727 
line, Zri A6143.180 line and Sii A6145.020 line. 



technical aspects of the spectrograph can be found in Zhao 
& Li (2001). 

The data reductions were carried out with a standard 
MIDAS package. The procedure is: order identification, 
background subtraction, flat-field correction, order extrac- 
tion, wavelength calibration, radial velocity shift correc- 
tion, and continuum normalization. Bias, dark current 
and scattered light corrections are included in the back- 
ground subtraction. The pixel-to-pixel sensitivity varia- 
tions were corrected by using the flat-field. The wave- 
length calibration was based on thorium-argon lamp spec- 
tra. In Fig. 1, we present a portion of spectrum of the 
sample star HD 98839 in the region near Ball A6141.727 
line. 

The measurements of spectral line EWs were done by 
applying two different methods: direct integration of the 
line profile and Gaussian fitting. The latter is preferable in 
the case of faint lines, but unsuitable for the strong lines 
in which the damping wings contribute significantly to the 
equivalent width. The final EWs are weighted averages of 
these two measurements, depending on the line intensity 
(see Zhao et al. 2000 for details). 

3. Stellar atmospheric parameters 

Model atmospheres are specified by four parameters: effec- 
tive temperature, surface gravity, overall metallicity and 
microturbulent velocity. 

Effective temperature T^g is determined from the uvby 
system index (6 — y) and [Fe/H] using the empirical cali- 
bration of Alonso et al. (1999). This calibration is based on 
a large sample of field and globular cluster stars roughly 
covering spectral type from FO to K5, and derived by ap- 
plying the infrared flux method. The photometric data of 
sample stars are taken from Olsen (1993) and Hauck & 
Mermilliod (1998). The error of the photometric data is 
a{b - ?/)=0.004 mag. Adopting the mean fT(Fe/H)=0.14 
dex from the spectroscopic analysis, the statistical error 
of Teff is estimated to about ±50 K. Considering a possi- 
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Table 1. Basic data of the saniple stars 



Star name 


Sp. 


V'mag 


B-V 


b-y 


TT 




BC 


M/Mq 












(mas) 








HD 26886 


G5 


7.00 


0.92 


0.596 


2.74 


1.05 


-0.335 


2.78(+0.75 -0.78) 


HD 27271 


G5 


7.52 


1.00 


0.591 


6.01 


1.13 


-0.306 


1.90(+0.25 -0.20) 


HD 50082 


GO 


7.43 


1.03 


0.611 


4.71 


0.99 


-0.367 


1.32(+0.28 -0.22) 


HD 98839 


Gsn 


4.99 


0.98 


0.610 


6.63 


0.63 


-0.308 


3.62(+0.45,-0.00) 



ble error of ±80 K in the cahbration, and other factors, we 
estimate the error in TcS as high as 100 K for our sample 
stars. 

The previous method to determine gravity logg re- 
quires that Fei and Fell Unes give the same iron abun- 
dance. But it is well known that the derivation of iron 
abundance from Pel and Fell lines may be affected by 
many factors such as unreliable oscillator strengths, pos- 
sible non local thermodynamic equilibrium (non-LTE) 
effects and uncertainties in the temperature structure of 
the model atmospheres. From Hipparcos parallaxes, we 
can determine more reliable gravities. Using the relations: 

log^ = log^ + Alog^^ + 0.4(Mboi - Aftoi,©) 

9q Jeff,© 

and 

Mboi = V -I- 5 + blog-K + BC, 

where M is the stellar mass, Mboi the absolute bolomet- 
ric magnitude, V the visual magnitude, BC the boloniet- 
ric correction, and tt the parallax. We adopt solar value 
log 50=4.44, Tcff,0=577O K, Mboi,0=4.75 mag. The par- 
allax TT is taken from the Hipparcos Satellite observations 
(ESA 1997). Stellar mass was determined from the po- 
sition of the star in Afboi^logT'off diagram. The corre- 
sponding errors were estimated by changing Tcff, [Fe/H] 
and A/boi according to its derivations (more details can be 
found in Chen et al. 2000). We adopt the stellar evolu- 
tion tracks given by Girardi et al. (2000). The bolometric 
correction, _BC, is determined using the empirical cali- 
bration of Alonso et al. (1999). The uncertainty of \ogg 
determined by this method mainly comes from the un- 
certainties in the stellar distance, stellar mass, effective 
temperature and bolometric correction. We estimate the 
error in \ogg to be in the range of 0.1—0.3 dex for our 
sample stars (Table 2 and the Col. 4 in Table 5). 

The initial nietallicities of the program stars were 
taken from the literature if available. Otherwise, we ob- 
tain an initial value by judging from the spectra and the 
color index. The final model metallicity is derived from 
consistency with the other parameters in the abundance 
calculation. 

Microturbulent velocity S^t was determined from the 
abundance analysis by requiring a zero slope of [Fe/H] vs. 
EWs. The values of the sample stars are given in Table 2. 
The error of is about 0.3 km s^^. 

Some characteristic parameters of four barium stars 
are presented in Table 1. The columns are: star name. 



Table 2. Atmospheric parameters of the sample stars 



Star name 


Tcff 


log 3 


it 


[Fe/H] 




(K) 




(km s"^) 




HD 26886 


4802 


2.22±0.3 


1.2 


-0.28±0.16 


HD 27271 


4874 


2.98±0.2 


1.4 


-0.06±0.14 


HD 50082 


4730 


2.50±0.2 


1.5 


-0.40±0.14 


HD 98839 


4866 


2.33±0.1 


1.6 


-f0.16±0.16 



spectral type, V magnitude, {B — V), (6 — y), parallax 
TT, uncertainty of parallax, bolometric correction and stel- 
lar mass. The atmospheric parameters of the sample stars 
are summarized in Table 2 including star name, effective 
temperature, surface gravity, microturbulent velocity and 
metallicity. 

4. Model atmospheres and spectral lines 

We adopt the model atmospheres generated by ATLAS9 
code (Kurucz 1993) to do the abundance analysis: 
these are LTE, plane-parallel, line-blanketed models. 
Abundances are determined by using the input atmo- 
spheric parameters given in Table 2 and the measured 
EWs. We usually choose the lines with EWs=20'-120 
mA, because the weaker lines would increase random er- 
ror and possibly some systematic overestimates, while 
the stronger lines are very sensitive to the microtur- 
bulence and damping. We use the updated oscillator 
strengths of spectrum lines taken from the NIST database 
(http://physics.nist.gov). But the log gf values of n- 
capture process elements are taken from the related ref- 
erences (Hannaford et al. 1982; Biemont et al. 1981; 
Weise & Martin 1980; Luck & Bond 1991; Biemont et 
al. 1982). All of the lines used in determining abundances 
are given in Table 3, including line and wavelength, exci- 
tation potential x, oscillator strengths logy/ and EWs. 
For Ban A5853.688 line in HD 98839, we find its EW is 
171 mA, which is much lower than 259 mA obtained by 
Fernandez- Villacanas et al. (1990). The EW values of Ban 
AA6141.727, 6496.908 lines in this star are higher than 200 
mA obtained in this study. 

5. Abundances and their uncertainties 

5.1. Abundance results 

The final abundances loge (in the usual scale logAr(H) 
= 12.0) and the corresponding induced [A'/Fe] values are 
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Table 3. Observed ion lines 



EWs (mA) 



AC 4) 


Ion 


Y CeVi 




HD 26886 


HD 27271 


HD 50082 


HD Q883Q 


UOVJU. 1 




4.61 


X . VJU 


7(1 

1 VJ 


77 






fiSriQ 994 


Fe I 


3.88 


— 1.84 


76 








UOUU.UC/U 


Fe I 


4.29 


X .00 


fi9 




fil 

VJ X 


7fi 

f VJ 




Fe I 


4 iif; 


— n fifi 

VJ.UVJ 






81 

X 






Fe I 


4 iif; 




QQ 








U U < J . u o u 


Fe I 


4 

4:.VJ^J 


— n 7^ 

VJ. 1 




8^ 
00 


7n 

J VJ 






Fe I 


2.45 


— 9 QQ 


Q4 




102 




t;q97 707 


Fe I 




— 1 flQ 


fi4 


fiQ 


47 


7fi 

f VJ 




Fe I 




— 1.41 




64 

VJ1 




8^1 
00 




Fe I 


4.65 


—0.23 


99 


112 


97 






Fe I 




— 1 1 7 

X . X ( 


Q7 










Fp I 


VJ.OVJ 


—4 fifl 

UVJ 


104 








fifin"^ 099 


Fe I 


3.88 


— 1.12 


106 


113 


104 






Fe I 


4 07 


— 1 flQ 


8fi 

OVJ 


Q4 


Q^ 

CJO 


120 


uvjuu.u ±o 


Fp I 
jr e I 


4 7*^ 


— n 4fi 

VJ.1VJ 


SI 

X 


Qn 

y VJ 


SS 
00 


114 


fin7Q nifi 


Fe I 


4.65 


— 1.12 




72 








Fc I 


4.61 


— 1 llfl 

X . UVJ 


■^fi 

OVJ 


■in 

OVJ 


41 


fi9 

VJZ 


fiflQfi fi71 


Fp I 
J? e I 


o. y o 


— 1 Q'^ 


00 


fiQ 
oy 




7Q 
f ^j 


fil fil fi9^ 

U -I- U -I- -U^O 


Fe I 


2.18 


— ^ 98 


81 

X 




Qn 

c/VJ 


1 1 fi 

X X VJ 


U -L U t J . O U O 


Fr> T 


4.14 


— 1 47 

X .4: J 


fill 

VJU 




fiQ 

VJc/ 




fil 7"^ "^41 


Fp I 
x" e I 


2.22 


— 9 SS 
^.00 


1 flQ 
xvjy 


1 1 Q 
X X y 


1 1 fi 
X X VJ 




fil sn 9nQ 


Fe I 


9 7^ 


— 9 fiS 
^. uo 


SS 

00 


1 1 ^ 

X xo 






fil S7 QQ'^ 


Fp I 


Q4 


— 1 79 




fil 

VJ X 


7n 

1 VJ 


Q'^ 
y 


fi9nn ^91 


Fp I 
e I 


9 fil 


— 9 44 




lis 

X xo 


1 14 
X x^ 




fi91 ii 14Q 


Fe I 


4.19 


— 1 1 ^ 

X . xo 


1 flQ 

XVJiJ 








fi99Q 9*^9 


Fc I 


2.84 


—2.81 


66 






108 


fi9^9 fi4S 


Fe I 


o.uu 


— 1.22 


112 










Fe I 


2.22 


— ^ 97 


Qfl 

C/VJ 








fi'i99 fiQ4 


Fe I 


9 iiQ 


—2.45 






1 1 <^ 

X X 




uoou.oou 


Fe I 


fiC) 


— n 86 

VJ.OVJ 


122 








UOOU. I uu 


Fe I 


4.19 


— 1.29 


77 




89 
0^ 


112 


fi4S1 87S 

U10 1 .0(0 


Fe I 


2.28 


— 9 Q7 


102 


1 1Q 

X X CI 


1 1 n 

X X VJ 




fifil S ^7^ 

Uu 1 0.O ( O 


Fe I 


9 8^ 


—2.45 




1 n^ 

XVJO 






U y i . t J U -L 


Fp I 
X e I 


4 sn 


— 1 fifi 

X . VJVJ 


00 






74 


UUUc/. 1 ±0 


Fe I 


9 iifi 


— 9 fifi 

^ .UVJ 


QQ 




1 1 ^ 

X X 




fi7n'^ "176 

U f UO.U f u 


Fp I 
X e I 


9 7fi 


— Ifi 

. X VJ 


74 


S'^ 

00 


74 


1 1 

X X 




Fe I 


4.61 


— 1.00 


60 










Fe I 


4.64 


— 1.58 


41 






60 


fi7f;n ifi4 


Fp I 
X e I 


2.42 


—9 fin 

^ . VJVJ 


104 




111 




fi7'i9 71fi 


Fe I 


4.64 


— 1.20 


50 








uouu.ouu 


Fe I 


9 7^ 


—3.21 


fifi 

VJVJ 


7^1 

f 


fiS 

VJO 


QQ 


fiSI fl 9fi7 


Fe I 


4.61 


—0.99 


62 


72 




89 


uo^o.uyu 


Fp I 
X e I 


4 fi4 


— Q9 

VJ • 


74 


SS 
00 




1 nfi 

XVJU 


UOOC/.OOU 


Fe I 


9 =lfi 








7n 

f VJ 


mil 

XVJO 


fiS49 fiSQ 


Fe I 


4 g4 


— 1.32 






44 




UO'iO.UUU 


Fe I 




— n Q^ 

VJ . C/O 


7Q 




78 
f 


1 1 ^ 

X X 


UOUO. 1 uu 


Fe I 


4.61 


— n Q^ 

VJ . C/O 




8^ 
00 






7071 Sfifi 
( VJ ( ± .ouu 


Fe I 


4.61 


— 1 7fl 

X . ( VJ 


"^fi 

OVJ 


47 


^7 
f 


00 


71 1 9 1 7n 


Fe I 


9 QQ 


—9 QQ 


61 

VJ X 




fiS 

VJO 




71 ^9 QS^i 


Fe I 


4 n7 


—1 fi^ 

X .VJO 


fin 

VJVJ 






QS 

iJO 


791 Q fiSfl 


Fe I 


4 07 


— 1 "^fl 
X .00 


fif; 

VJU 


71 

f X 




Qn 

iJVJ 


7^nfi ii7fl 

f OUU.U I VJ 


Fe I 


4.18 


— 1 74 

X . ( 1 


fi7 
f 








741 S fi79 


Fe I 


4.14 


—1.38 


67 




73 


108 


f UOO. 1 <yU 


Fe I 


^ 09 


— 1 88 
X .00 


120 








7723.210 


l^e I 


2.28 


—3.62 


73 






112 


8365.640 


Fei 


3.25 


-2.04 






115 




8515.122 


Fei 


3.02 


-2.07 


111 








8598.836 


Fei 


4.39 


-1.09 




90 






8699.461 


Fei 


4.95 


-0.38 


75 


89 




100 


5991.378 


Fell 


3.15 


-3.56 




52 


47 


79 


6149.249 


Fell 


3.89 


-2.72 


52 




50 
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Table 3. Continued 



EWs (mA) 



AC 4) 


Ion 


Y feVl 
}L V^^ ) 




HD 26886 


HD 27271 


HD 50082 


HD 9883Q 


fi41 fi Q9S 


Fe II 




— 9 74 


4Q 


OO 


47 


7S 




Fe II 


2.89 


—3.58 




55 


61 






Fe II 




— 9 07 

Z .U ( 


Qii 

C/O 


78 
f o 






UU 1 u.uoo 


Fe II 


9 8Q 


— 3 1 fl 

O. J-U 


88 

OO 






1 in 

X xu 


771 1 7^1 


Fp ti 

jr c 11 


^ Qfl 


— 9 47 




OO 






7771 


O T 


9.14 


fl 33 


ii9 
oz 


111 

O 1 




fi4 


7774 1 77 


D T 


9.14 


fl 1 Q 

U. J. c/ 




llfi 

ou 




fif; 

UO 


fil fi4 9^n 


Nai 


2.10 


— 1 ii7 


iiO 
ou 


so 

ou 


ii7 

O f 


in7 

xu f 


filfifl 7^1^ 

U lUU. f OO 


Nai 


9 1 n 


— 1 93 


79 








S71 7 8'^'^ 


Mo' T 


9.53 


—0.97 


77 


97 


40 


100 




All 


^ 14 


— 1 87 

1 .0 ( 


38 

OO 




37 

O f 




7S^ii ^1 7 


All 


4.02 


— fl iiS 
u.uo 


4fi 

'iU 


70 
f u 




S3 

OO 


1 OOU . -L O 


All 


4.02 


—0.40 


62 


85 








Si T 


o . uo 


— 1 1 Q 


Q4 


1 1 n 




124 


fil n9fi 


Si T 


^ fil 




40 

'iU 


ii7 

O 1 


OO 




fil 49 4Q4 


Si I 


5.62 


— 1.48 


43 








fi14'i 09(1 


Si T 


fi9 
u . uz 


— 1 4"^ 








7S 


74n'i 7Qn 


Si I 


5.61 


—0.68 


100 


90 




139 


741 ^ QfiS 


Si T 


fil 


— fl 71 

U. f 1- 






111 




7sno nnn 

1 ouu.uuu 


Si T 


fi 1 8 


— n 78 

U. 1 o 


iiS 

OO 








7Q'^9 "^iil 


Si I 


5.96 


—0.35 


95 






116 


fiSfi7 ii79 


Ca I 


9 


— 1 fil 

1 .U 1 






3fi 
ou 




fil fil 9Q"i 


Cai 


2.52 


— 1.19 


98 








filfifi 44n 


Cai 


2.52 


— 1 1 Q 


Q4 


1fl4 


1fl4 


120 


filfiQ n44 


Ca I 


2.52 


— fl Sfl 

u.ou 


1 


132 




1 flfi 

X ou 


fi4'i'i finii 


Ca T 


2.52 


— 1.29 




Q7 


8*^ 

OO 




fi471 fifiS 


Ca I 


2.52 


— n fiQ 


121 


1 3Q 


1 97 

X Z ( 


1 fin 

X ou 


fi94T fi9n 


Sr TT 

kJV- 11 


1 t1 


— n Qs 


77 








fifin4 fifin 


Sr IT 


1 ^fi 


— 1 1 fi 
J. . J.U 


71 




S3 

OO 


Q7 


fiSfifi 4fi1 


Ti I 


1 n7 


— n 84 


Q1 






132 




Ti I 


1 SQ 


— n 90 




S3 

OO 






fil 9fi 994 
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Table 3. Continued 



EWs (mA) 

A(A) Ion X (eV) logg/ Rcf." HD26886 HD27271 HD50082 HD98839 



6435.000 


Yi 


0.07 


-0.82 


h 


25 


47 


51 




6127.460 


Zrl 


0.15 


-1.06 


b 




44 


49 


— 


6134.570 


Zrl 


0.00 


-1.28 


b 




36 


44 




6140.460 


Zrl 


0.52 


-1.41 


b 


5 


16 


17 




6143.180 


Zrl 


0.07 


-1.10 


b 


27 


44 


69 




5853.688 


Ban 


0.60 


-1.01 


wm 




158 


200 


171 


6141.727 


Ban 


0.70 


-0.08 


wm 


255 


239 


370 


246 


6496.908 


Ban 


0.60 


-0.38 


wm 


245 


235 


341 


257 


6390.480 


Lan 


0.32 


-1.45 


lb 


48 


48 


82 


51 


6645.110 


Eun 


1.37 


0.20 


b2 


30 


28 


41 


45 



References — h: Hannaford et al. (1982); b: Biemont et al. (1981); wm: 
Weise & Martin (1980); lb: Luck & Bond (1991); b2: Biemont et al. (1982). 



presented in Table 4, where N refers to the hne number 
used in ealeiihition. The solar abundanees are taken from 
Grevesse & Sauval (1998). 

Fig. 2(a)- (d) shows the detailed abundance results of 
the sample stars. The corresponding errors are taken from 
Col. 6, atot, in Table 5. Detailed analysis about errors is 
discussed in next section (Sect. 5.2). 

5.2. Errors in resulting abundances 

Uncertainties of the abundances mainly come from equiv- 
alent widths measurement, log gf values and stellar atmo- 
spheric parameters. 

Uncertainties in the EWs are set essentially by the S/N 
and the resolution of the spectra. An expression for the 
intrinsic accTiracy of an equivalent width is provided by 
Cayrel (1988). Applied to our spectra, having R=40,000 
and typical S/N of 150, we expect micertainties in the 
EWs of about 2 mA. We also consider the errors obtained 
by measuring EWs of spectral lines several times in the 
total uncertainties. For an element represented by N lines, 
the error is decreased by a factor \/7V. Thus the errors in 
element abundances from EWs are calculated and illus- 
trated on Col. 2 in Table 5, asw/VN- 

The errors caused by atmospheric parameters can be 
obtained through changing the parameters Tcfr, log g, 
individually. Column 3, 4, 5 in Table 5 show the effects on 
the derived abundances by changing -1-100 K in TeS, +0.3, 
+0.2 or +0.1 dex in logg, +0.3 km s~^ in 

Model metallicity uncertainties scarcely affect the ele- 
ment abundances, and they can be ignored. 

The uncertainty in (//-values of the Fe I and Fe ii lines 
is about 0.05 dex. For the elements up to Ni, the uncertain- 
ties are about 0.02—0.10 dex, and for the heavier elements, 
the errors are ~0.20 dex (Smith et al. 1995; Lambert et 
al. 1996; Junqueira & Pereira 2001). 

The combined abundance errors caused by the above- 
mentioned uncertainties in (//-values, EWs-values and 
stellar atmospheric parameters, are the total errors in 
abundances, which are given in the last columns in Table 
5.(l)-(4) hy atot. 



5.3. Abundance pattern analysis 

Among the four sample stars, HD 98839 shows the high- 
est metallicity with [Fe/H]=+0.16, which is slightly higher 
than solar metallicity. HD 50082 shows the lowest metal- 
licity with [Fe/H]= 0.40. The other two stars ([Fe/H] is 
-0.06 for HD 27271 and -0.28 for HD26886) are nearer to 
solar metalUcity than HD 98839 and HD 50082. However, 
it is clear that all of them are disk stars. 

The [q/Fc] of sample stars are similar to solar. And 
it is clear that [a/Fe] patterns are lower than those of 
the two halo Ba stars ([Fe/H] =-1.72, -1.43) studied by 
Jimqueira & Pereira (2001). Their two halo Ba stars show 
the overabundant [a/Fe] relative to solar. The reason is 
that massive stars produce these a elements, which have 
been ejected into interstellar medium with the Type II 
supernova explosion of the massive stars since the early 
stage of our Galaxy. Most of element Fe is produced 
through Type la supernova explosion, which occurs in 
the close binary system composed by the longer lifetime, 
intermediate- and low-mass stars (Timmes et al. 1995; 
Samland 1998; Liang et al. 2001). Thus, there is abim- 
dant element Fe in our Galaxy at its present stage. The 
[a/Fe] patterns of our sample Ba stars also show that all 
of these four stars are disk stars. In addition, it should 
be noticed that the [0/Fe] values of our three samples 
will be increased by 0.14 when we use the recent solar 
O abundance loge(0)=8.69 obtained by AUende Prieto 
et al. (2001) based on [Oi] A6300 forbidden line. That is, 
the [0/Fe] ratio in HD 26886, HD 27271 and HD 98839 are 
0.50, 0.53 and 0.29, respectively. But this increase value is 
still smaller than the corresponding abundance errors in 
the three sample stars: 0.22, 0.18 and 0.21, respectively. 

For Na and Al, these two odd elements, their abun- 
dances are similar to Mg abundance, because these nuclei 
are thought to be synthesized in the same process, namely, 
carbon burning in massive stars. 

The iron group elements follow Fe closely though with 
some scatter. Generally, Cr and Mn are slightly overabun- 
dant, and Ti, V are slightly underabundant. Especially, 
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Fig. 2. (a) (the top-left panel) Abundance pattern of HD26886 with the error bars taken from Col.6 in Table 5(1); 
(b) (the top-right panel) Abundance pattern of HD27271 with the error bars taken from Col.6 in Table 5(2); (c) 
(the bottom-left panel) Abundance pattern of HD50082 with the error bars taken from Col.6 in Table 5(3); (d) (the 
bottom-right panel) Abundance pattern of HD98839 with the error bars taken from Col.6 in Table 5(4). 



Ni, the closest element to Fe, is very similar to the Fe 
abundance in the sample stars. 

However, it is not clear whether the larger scatter of 
a elements and iron elements in HD 50082 is the natural 
results of its strong Ba property or not. 

The most important results are the obvious overabun- 
dance of s-proccss elements {Z >56), Y, Zr and Ba. La, 
Eu (two peaks), in all of the four Ba stars (without appro- 
priate Y and Zr lines are used in HD 98839 due to spectral 
qualities). Especially, the second peak s-process elements, 
Ba and La, are in the highest abundances. The origin of 
the two peaks is the result of neutron magic number 50 
and 82 nuclei on the path of s-process nucleosynthesis oc- 
curred in the interiors of AGB stars (Clayton et al. 1961; 
Seeger et al. 1965), the progenitor companions of these 
Ba stars. There may be a third peak on clement Pb corre- 
sponding to the neutron magic number 126. But it needs 
high quality blue spectra to get the Pb i AA3683, 4057 lines 
to obtain the Pb abundance. 

To further degree, in the three sample Ba stars, the 
first peak s-elements, Y and Zr, show lower abundances 



than the second peak s-elements Ba, La etc. This trend 
also was shown in some other Ba stars given by Zacs 
(1994). This [hs/ls]>0 trend ("hs" refers to the second 
peak elements, "Is" refers to the first peak elements) shows 
that the neutron exposures experienced in their AGB pro- 
genitor companions are strong enough to produce the 
more abundant "hs" elements than the "Is" elements. If we 
plot our Ba samples and the corresponding AGB progeni- 
tor companions on Fig. 1 of Liang et al. 2000), or Fig. 6 of 
Busso et al. (1995) ([Is/Fe] vs. [hs/ls]), these stars should 
be in the [hs/ls] >0 regions in the two figures. Our theo- 
retical calculations given in next section certify this point. 

This discrepancy between the two peak elements may 
be caused by the different neutron exposures in the com- 
panion AGB stars. Namely, the higher neutron exposure 
in the interior of AGB stars benefits to the production of 
the heavier s-process elements (Liang et al. 2000; Busso et 
al. 1995). 

Comparing the abundance patterns of the two "typi- 
cal" Ba star HD 27271 and HD 50082, we notice that the 
differences between the two peak elements (Y, Zr and Ba, 
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Table 4. Element Abundances in the usual scale logX(H) = 12.0 and in the notation [X/Fe] in the four Ba stars 



[Fe/H] 
Ion 


N 


HD 26886 
-0.28 
loge 


[X/Fe] 


N 


HD 27271 
-0.06 

loge 


[X/Fe] 


N 


HD 50082 
-0.40 
loge 


[X/Fe] 


N 


HD 98839 
+0.16 
loge 


[X/Fe] 


Fel 


46 


7.23 




26 


7.45 




29 


7.11 




27 


7.67 




Fell 


5 


7.35 


— 


6 


7.49 


— 


5 


7.26 


— 


4 


7.82 


— 


Oi 


1 


8.91 


+0.36 


2 


9.16 


+0.39 




— 


— 


2 


9.14 


+0.15 


Nai 


2 


5.97 


-0.08 


1 


6.43 


+0.16 


1 


6.01 


+0.08 


1 


6.78 


+0.29 


Mgl 


1 


7.48 


+0.18 


1 


7.62 


+0.10 


1 


6.92 


-0.26 


1 


7.72 


-0.02 


All 


3 


6.21 


+0.02 


2 


6.44 


+0.03 


1 


6.22 


+0.15 


1 


6.61 


-0.02 


Si I 


6 


7.49 


+0.22 


3 


7.70 


+0.21 


2 


7.69 


+0.54 


5 


7.90 


+0.19 


Cai 


4 


6.18 


+0.10 


4 


6.29 


-0.01 


6 


5.97 


+0.01 


3 


6.53 


+0.01 


Sen 


2 


2.82 


-0.07 


1 


3.11 


0.00 


1 


2.93 


+0.16 


1 


3.23 


-0.10 


Til 


6 


4.47 


-0.27 


4 


4.79 


-0.17 


4 


4.48 


-0.14 


3 


4.95 


-0.23 


Vi 


1 


3.46 


-0.26 


2 


3.61 


-0.33 


2 


3.22 


-0.38 


1 


3.92 


-0.24 


Cri 


1 


5.52 


+0.13 


2 


5.87 


+0.26 


2 


5.43 


+0.16 


2 


6.06 


+0.23 


Mm 


2 


5.28 


+0.17 


1 


5.35 


+0.02 


1 


5.22 


+0.23 


1 


5.77 


+0.22 


Nil 


14 


5.94 


-0.03 


14 


6.18 


-0.01 


5 


5.75 


-0.10 


9 


6.37 


-0.04 


Yl 


1 


1.91 


-0.05 


1 


2.65 


+0.47 


1 


2.30 


+0.46 








Zrl 


2 


2.42 


+0.10 


4 


2.95 


+0.41 


4 


2.80 


+0.60 








Ban 


2 


2.59 


+0.74 


3 


2.74 


+0.67 


3 


2.87 


+1.14 


3 


2.85 


+0.56 


Lan 


1 


1.56 


+0.67 


1 


1.92 


+0.81 


1 


2.12 


+1.35 


1 


1.73 


+0.40 


Eun 


1 


0.49 


+0.26 


1 


0.81 


+0.36 


1 


0.74 


+0.63 


1 


0.94 


+0.27 



La) are higher in HD 50082 than in HD 27271. It may be 
that the companion AGB star of HD 50082 experiences 
stronger neutron exposure in s-process nucleosynthesis. As 
a result, the abundances of all of the s-process elements 
of HD 50082 are higher than those of HD 27271. Our the- 
oretical calculations certify this analysis (see Sect. 6). 

As for calling HD 50082 and HD 27271 as "typical" Ba 
stars, the main reason is their masses. According to the 
scale defined by Warner (1965), there are three mild Ba 
stars in our samples: HD 26886, HD 27271 and HD 98839, 
and one strong Ba star HD 50082. The average masses 
are 1.9M0 and 1.5M0 respectively for mild Ba stars and 
strong Ba stars (Jorissen ct al. 1998). We obtain masses 
of 1.9OM0 and 1.32M0 for HD 27271 and HD 50082 re- 
spectively, which are consistent with the mass judgments 
well. 

HD 26886 and HD 98839 also show over-abundant Ba 
and La, but the abundances are lower than those of 
HD 50082 and HD 27271. The elements Y, Zr of HD 26886 
are weakly over-abundant. These results can be under- 
stood basically by the weaker overabundances of heavy el- 
ements in their AGB progenitor companions, by the lower 
neutron exposures occurred in the AGB progenitors. 



6. Comparison with wind accretion model results 

What is the reason to cause the heavy-element overabun- 
dances of Ba stars? It is generally believed that the over- 
abundances are caused by binary accretion, namely, Ba 
stars accret the ejected materials from their compan- 
ions, the former AGB stars (the present white dwarfs), 
which synthesized these heavy elements by themselves and 



ejected the elements into interstellar medium through stel- 
lar wind (Liang et al. 2000 and references therein). 

In our previous paper (Liang et al. 2000), we calcu- 
lated the s-process nucleosynthesis of AGB star with 3M0 
and solar metallicity. At the same time, we set up an an- 
gular momentum conservation model of wind accretion 
in binary systems. Using such model we calculated the 
heavy-element abundances of Ba stars, and successfully 
explained the observed abundances of some samples taken 
from Zacs (1994) and Busso et al. (1995). 

We try to explain the abundance patterns of our sam- 
ple stars using our theoretical method. The calculation of 
the heavy clement abundances of Ba stars is made in two 
separate steps. Firstly, adopting the theory of s-process 
nucleosynthesis and the latest TP-AGB model (Straniero 
ct al. 1995; Straniero ct al. 1997; Gallino ct al. 1998; Busso 
et al. 1999), we calculate the overabundances of the intrin- 
sic AGB star at each ejection. Then, combining the accret- 
ing matter predicted by the model of wind accretion on 
successive occasions and mixing, we calculate the heavy- 
element overabundances of the barium star. The s-process 
nucleosynthesis scenario and the change equations of or- 
bital semi-major axis A and eccentricity e, and j^^, 
were given in Liang et al. (2000) and Liu et al. (2000). 

The standard case of wind accretion is: Mi,o=3.OM0, 
M2,o=l-3M0, Ve3=l5 km s~^ {Mi^ is the main sequence 
mass of the intrinsic AGB star, the present white dwarf, 
in the binary system; Af2,o is the corresponding mass of 
the present Ba star; Vej is the wind velocity). The used 
accretion rate is 0.15 times of the Bondi-Hoyle's accretion 
rate (Liang et al. 2000; Boffin & Zacs 1994). 

We try to explain the observed abundances of our three 
Ba stars, HD 50082, HD 27271 and HD 26886, from theo- 
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HD 26886 





Vn 


AT <T 

^-t eff 


Alogg 




^tot 




+100 K 


+0.3 


1 n o 1 ... — 1 

+0.3 Km s 




Fei 


0.02 


0.08 


0.01 


—0.14 


0.16 


Fe II 


0.07 


—0.08 


0.14 


—0.14 


0.22 


O I 


0.12 


—0.15 


0.11 


—0.03 


0.22 


Nai 


0.07 


0.08 


—0.01 


—0.04 


0.11 


Mg I 


0.10 


0.03 


—0.03 


—0.04 


0.11 


All 


0.05 


0.08 


—0.01 


—0.03 


0.10 


Si I 


0.05 


—0.01 


0.04 


—0.06 


0.09 


Cai 


0.10 


0.11 


—0.03 


—0.16 


0.22 


Sen 


0.10 


—0.01 


0.12 


—0.11 


0.19 


Ti I 


0.05 


0.15 


—0.01 


—0.10 


0.19 


V I 


O.lo 


n 1 c 
O.lo 


0.00 


—0.09 


U.zz 


Cn 


0.19 


0.12 


-0.02 


-0.18 


0.29 


Mm 


0.15 


0.12 


-0.04 


-0.20 


0.28 


Nil 


0.04 


0.06 


0.03 


-0.14 


0.16 


Yi 


0.05 


0.24 


0.08 


-0.02 


0.27 


Zrl 


0.04 


0.19 


0.03 


-0.01 


0.21 


Ban 


0.08 


0.03 


0.01 


-0.07 


0.14 


Laii 


0.11 


0.02 


0.13 


-0.09 


0.21 


Euii 


0.08 


-0.01 


0.13 


-0.05 


0.18 


Table 5.(3). 


Abundances uncertainties of HD 50082 




HD 50082 


Ion 






Alog (/ 


Mf 






+ 100 K 


+0.2 


+0.3 km s"' 




Fei 


0.03 


0.07 


0.01 


-0.12 


0.14 


Fell 


0.05 


-0.09 


0.10 


-0.08 


0.16 


Nai 


0.08 


0.08 


-0.01 


-0.03 


0.12 


Mgi 


0.06 


0.02 


0.00 


-0.01 


0.07 


All 


0.06 


0.07 


-0.00 


-0.02 


0.10 


Si I 


0.09 


-0.02 


0.03 


-0.06 


0.11 


Cai 


0.06 


0.10 


-0.01 


-0.10 


0.15 


Sen 


0.15 


-0.01 


0.08 


-0.10 


0.19 


Til 


0.06 


0.16 


0.00 


-0.07 


0.18 


Vi 


0.09 


0.16 


0.01 


-0.09 


0.20 


Cn 


0.10 


0.10 


-0.01 


-0.12 


0.18 


Mm 


0.19 


0.11 


-0.04 


-0.18 


0.28 


Nil 


0.05 


0.06 


0.04 


-0.08 


0.12 


Yi 


0.07 


0.20 


0.06 


-0.04 


0.23 


Zrl 


0.04 


0.18 


0.03 


-0.03 


0.21 


Ban 


0.08 


0.03 


0.00 


-0.07 


0.14 


Laii 


0.16 


0.03 


0.08 


-0.16 


0.26 


Euii 


0.08 


-0.01 


0.09 


-0.05 


0.15 



retical calculations. Notice that there is no eccentricity of 
HD 98839 was observed, so we do not calculate the the- 
oretical abundances to fit its observations. In our calcu- 
lations, we try to make the calculated orbital period P 
and eccentricity e match the observations of the samples. 
Tabic 6 lists the observed orbital elements and their Ba 
classes of the samples taken from Jorissen et al. (1998), 
and their masses obtained by us (also in Table 1). 

As the studies show, the Ba stars with orbital pe- 
riod P >1600 days can be formed through wind accretion 
(Zhang et al. 1999; Liang et al. 2000; Liu et al. 2000). 
Jorissen et al. (1998) suggested the corresponding period 
is 1500 days. Possibly, the Ba stars with lower orbital pe- 



HP 27271 

1^ ATefi Alogs Ai; ^ 



+100 K +0.2 +0.3 km g-^ 



Fei 


0.03 


0.05 


0.01 


-0.13 


0.14 


Fell 


0.07 


-0.03 


0.17 


-0.03 


0.19 


Oi 


0.08 


-0.14 


0.07 


-0.03 


0.18 


Nai 


0.11 


0.08 


-0.02 


-0.05 


0.15 


Mgi 


0.09 


0.03 


-0.03 


-0.03 


0.11 


Ah 


0.07 


0.05 


-0.02 


-0.04 


0.11 


Si I 


0.07 


-0.02 


0.02 


-0.07 


0.10 


Cai 


0.09 


0.11 


-0.03 


-0.14 


0.20 


Sen 


0.13 


-0.01 


0.08 


-0.11 


0.19 


Til 


0.06 


0.14 


-0.00 


-0.08 


0.17 


Vi 


0.10 


0.16 


0.00 


-0.12 


0.22 


Cn 


0.11 


0.10 


-0.02 


-0.13 


0.20 


Mm 


0.18 


0.12 


-0.04 


-0.16 


0.27 


Nil 


0.04 


0.04 


-0.02 


-0.11 


0.12 


Yi 


0.07 


0.21 


0.02 


-0.04 


0.24 


Zn 


0.04 


0.19 


0.01 


-0.03 


0.21 


Ban 


0.09 


0.03 


0.01 


-0.11 


0.16 


La II 


0.10 


0.02 


0.09 


-0.07 


0.17 


Euii 


0.07 


-0.01 


0.09 


-0.04 


0.13 


Table 5.(4). 


Abundances uncertainties of HD 98839 




HD 98839 


Ion 




+ 100 K 


Alog q 
+0.1 


AO 
+0.3 km s"' 


ntof 


Fei 


0.03 


0.06 


0.01 


-0.14 


0.15 


Fell 


0.09 


-0.09 


0.05 


-0.15 


0.20 


Oi 


0.13 


-0.16 


0.04 


-0.05 


0.21 


Nai 


0.15 


0.10 


-0.11 


-0.08 


0.22 


Mgi 


0.11 


0.03 


-0.01 


-0.05 


0.13 


All 


0.11 


0.05 


-0.01 


-0.05 


0.13 


Si I 


0.07 


-0.02 


0.01 


-0.09 


0.11 


Cai 


0.13 


0.10 


-0.01 


-0.18 


0.24 


Sen 


0.16 


-0.01 


0.04 


-0.12 


0.20 


Til 


0.09 


0.16 


0.00 


-0.11 


0.21 


Vi 


0.10 


0.16 


0.00 


-0.13 


0.23 


Cn 


0.14 


0.11 


-0.00 


-0.16 


0.24 


Mm 


0.21 


0.11 


-0.01 


-0.22 


0.32 


Nil 


0.06 


0.04 


0.02 


0.01 


0.09 


Ban 


0.08 


0.03 


0.01 


-0.15 


0.17 


La II 


0.09 


0.02 


0.04 


-0.06 


0.12 


Eu II 


0.08 


-0.01 


0.04 


-0.06 


0.15 



riod form through other scenarios: dynamically stable late 
case C mass transfer or common envelope ejection. 

For our sample of Ba stars, the orbital periods of 
HD 50082 and HD 27271 are 2896 and 1693.8 days (>1600 
days) respectively. Their masses are consistent with the 
average masses of typical strong Ba stars (I.SMq) and 
mild Ba stars (I.9M0) (see Sect. 5.3), so we calculate the 
theoretical abundances of Ba stars using the standard case 
of our wind accretion model to compare with the observed 
patterns. Fig. 3(a), (b) show that the calculations (the 
solid lines) can fit observations very well within the errors 
for HD 50082 and HD 27271. The variable "a" represents 
the times of the corresponding standard exposures in the 
■^^C profile suggested by Gallino et al. (1998) (details can 
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Table 6. Masses and orbital elements of the sample stars 



Star name 


M/Mq 


P 


e 


Ba 






(days) 




classes 


HD 26886 


2.78(+0.75,-0.78) 


1263.2 


0.39 


mild 


HD 27271 


1.90(+0.25,-0.20) 


1693.8 


0.22 


mild 


HD 50082 


1.32(+0.28,-0.22) 


2896.0 


0.19 


strong 


HD 98839 


3.62(+0.45,-0.00) 


>11000 




mild 



be found in Liang et al. 2000). In summary, the higher 
a value reflects the higher neutron exposure occurred in 
interiors of AGB progenitor companions. 

Fig. 3(a), (b) show that the required neutron exposure 
characteristic parameter a is higher for HD 50082, a=2.6; 
while for HD 27271, a=1.6 is adequate. This difference re- 
flects the "strong" and "mild" Ba properties of each star. 

For HD 26886, since its mass is 2.78M0 (though with 
large error), we adopt 2.5M0 as the initial mass of a 
Ba star for a binary system in our calculation, instead 
of 1.3M0 in the standard case. However, the fit failed 
(Fig. 3(c), the bottom panel). Maybe this is due to the 
lower orbital period of HD 26886 (1263.2 days) than the 
lower limit of wind accretion scenario (1600 or 1500 days). 
Perhaps dynamically stable late case C mass transfer or 
common envelope ejection scenario is the formation source 
of HD 26886. However, more studies are needed to under- 
stand this suggestion. 

7. Discussions 

7.1. Stellar masses and mild Ba stars, strong Ba stars 

Adopting stellar evolution tracks given by Girardi et al. 
(2000) and some related parameters (see Sect. 3), we cal- 
culate the masses of the sample Ba stars (see Table. 1, 
6). Their masses are consistent with the statistical and 
theoretical results. Table 9 of Jorissen et al. (1998) shows 
that the average mass of mild Ba stars is 1.9Mq with 
O.6OM0 companion white dwarfs. For HD 27271, a mild 
Ba star, its mass is 1.9OM0, which fits the average I.9M0 
of mild Ba stars very well. As for HD 26886, the calculated 
2.78M0 is consistent with the average 2.3M0 of mild Ba 
stars with O.67M0 companion white dwarfs (Jorissen et 
al. 1998) within the errors. HD 98839 is discussed in next 
section. 

For the strong Ba star HD 50082, we obtain a mass of 
1.32 M0, which is consistent with the average mass I.5M0 
of strong Ba stars with O.6OM0 companion white dwarfs 
given by Jorissen et al. (1998). 

Besides the; different masses between the strong and 
mild Ba stars, their abundance patterns are different too 
(Sect. 5.3). What causes the discrepancy of abundances in 
these two Ba groups? Jorissen et al. (1998) carefully dis- 
cuss the possible reasons. Considering our results, we em- 
phasize and analyze two possible reasons again here: (i) a 
lower mass for the barium star results in a smaller dilution 
of the accreting materials in envelope of the Ba star; (ii) 
the lower metallicity results in larger heavy-element over- 




30 40 50 60 70 

Chorqe Z 




Charge Z 



HD 26886 




Charge Z 

Fig. 3. The fittings of the predicted to observed heavy 
element abundances of barium stars using our wind ac- 
cretion models, (a) for HD 50082 (the top panel); (b) for 
HD 27271 (the middle panel); (c) for HD 26886 (the bot- 
tom panel). 



abundances in the AGB progenitor. Our results are consis- 
tent with these two points. The strong Ba star HD 50082 
has the lowest mass and lowest metallicity in the four sam- 
ples. Thus it shows the highest heavy-element overabun- 
dances. 
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7.2. Inspiration from HD 98839: a "false Ba" star? 

For HD 98839, the atmosphere parameters, [Fe/H] and 
over-rich Ba abundance obtained by us are similar to 
those of other researchers (Pilachowski 1977; Fernandez- 
Villacahas et al. 1990; Mishenina 1996). We estimate its 
mass to be 3.62M0. Its high mass and heavy-element over- 
abundance pose some questions on Ba stars formation 
through binary accretion and also on nucleosynthesis sce- 
nario of AGB stars. 

The mass of HD 98839 is 3.62M0, higher than the aver- 
age masses of mild Ba stars (1.9 or 2.3M0, sec Sect. 7.1). 
And there is also a discrepancy between this value and 
the suggestion of Han et al. (1995), in which the authors 
assigned masses of less than 3M0 for all Ba stars. 

How can these results and discrepancies be under- 
stood? On one hand, this mild Ba star may have been 
classified erroneously as Ba star, namely it is "false Ba" 
star, and the real story is that it produced the over- 
abundant heavy elements by itself. However, according 
to the derived parameters, Mboi=— 1.2 and Teff=4866K, 
HD 98839 should not evolve to AGB stage. Generally, 
AGB stars have Mboi from —1.7 to —5.5 and Tes from 
3200 K to 3850 K respectively (Smith & Lambert 1985, 
1986, 1990). In this case, if the heavy element overabun- 
dances of HD 98839 is an intrinsic nucleosynthesis event 
within itself, the over-rich s-process elements is usually 
associated with mixing process at the helium core flare, 
which is consistent with the alternative origin for the 
Ban star phenomenon suggested by Malaney (1987) and 
Malaney & Lambert (1988). Also, Mennessier et al. (1997) 
suggested that some of their mild Ba stars with mass 
higher than 3M0 may be lower luminosity, post He-flare 
stars. 

Moreover, HD 98839 displays very long orbital period 
(>11000 days) up to the upper detected limit (Jorissen 
et al. 1998); and there is no observed binary orbital ec- 
centricity. Thus, it is not unrealistic to consider it as an 
independent star. 

However, if these high mass mild Ba stars are re- 
ally "true Ba" stars, and the binary accretions cause 
their heavy clement overabundances, it will challenge the 
present nucleosynthesis scenario of AGB stars. At present, 
the popular viewpoint is that 1— 3M0 low mass AGB stars 
with ^^C neutron source are the main nucleosynthesis sites 
of s-process elements, while the intermediate mass AGB 
stars with ^^Ne neutron source cannot fit to the obser- 
vations (Malaney 1987; Busso et al. 1995; Busso et al. 
1999: Liang ct al. 2000). However, if masses of Ba stars are 
higher than 3M0, their companions should be intermedi- 
ate mass stars on the main-sequence. How the intermedi- 
ate mass AGB stars result in the heavy element overabun- 
dance of these Ba stars? If this is the case, what kind of 
neutron exposure scenario works? Single exposure? More 
studies are needed on s-process nucleosynthesis scenario 
of AGB stars. Certainly, we need also the large sample of 
high mass mild Ba stars to check its abundance patterns. 



In addition to the high mass (3.62M0) of mild Ba star, 
HD 98839, obtained by us, there are also some other high 
mass samples. Smith ct al. (1980) calculated mass of C Cap 
is about 3.5M0, but with large error (-|-3.6,-1.8M0) (also 
see Boffin & Jorissen 1988). Mennessier et al. (1997) stud- 
ied higher mass Ba stars. They estimated that the masses 
are up to 7M0 for their group S and 4.5M0 for group C in 
which there are mild Ba stars. Monte Carlo simulations of 
Karakas et al. (2000) predict that the masses of Ba stars 
may be up to 6M0. 

Therefore, it is necessary to pay more attention on 
these high mass mild Ba stars. 

8. Summary 

We obtain abundances of a, iron peak and s-process el- 
ements in four Ba stars, HD 50082, HD 27271, HD 26886 
and HD 98839, using high resolution and high S/N spectra. 
The results show that all of samples are disk stars. They 
show the similar metallicities and the similar abundances 
of a-elments and iron peak elements to the solar values. 
All of them show s-process element overabundances rel- 
ative the solar. We obtain their masses using the highly 
precise Hipparcos data, precise photometric parameters 
and stellar evolution tracks. The mass results fit well the 
average characteristic masses of strong and mild Ba stars. 
And the high mass of mild Ba star HD 98839 pose more 
thoughts on Ba star phenomenon. 

Using our wind accretion model of binary systems and 
AGB stars s-element nucleosynthesis calculation, we con- 
clude that the wind accretion formation scenario is pos- 
sible for HD 50082 and HD 27271. Wind accretion cannot 
explain the observations of HD 26886 with the short or- 
bital period (P <1600 days). The mild Ba star HD 98839 
with high mass and very long orbital period may be ei- 
ther a "false Ba" star with the heavy elements enriched 
by itself or a "true Ba" star. 

Fortunately, though we only have four samples, they 
represent wide Ba star groups. From stellar masses, there 
is one strong Ba star among them, and also mild Ba stars; 
there are the members with typical mass of Ba stars, also 
there are the high mass mild Ba stars, which are contro- 
versial. From orbital elements, there are samples in typical 
wind accretion orbital period range, there is one sample 
with very long orbital period, up to the detected limit, and 
there is one sample with shorter value, which may corre- 
spond to other formation scenarios. So these samples are 
very interesting in understanding the abundance patterns 
and formation scenarios of Ba stars. 

Certainly, we need larger samples to understand the 
nature of Ba stars. At present, there is a large sample of Ba 
stars with orbital elements (Carquillat et al. 1998; Udry 
et al. 1998a, 1998b; Jorissen et al. 1998), but without the 
corresponding element abundances. It is very interesting 
and necessary to obtain their element abundances to be 
combined with the orbital elements. It is specially impor- 
tant to compare the related results of the two groups with 
different orbital periods: P <1600 and P >1600 days, and 
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to study the stars with very long orbital periods. These re- 
sults will help us to understand the various formation sce- 
narios of Ba stars. At the same time, it is important to ob- 
tain their masses, particularly, for the possible high mass 
mild Ba stars. Maybe these high mass mild Ba stars can 
give more clues on AGB stars nucleosynthesis. Aforcovcr, 
it will be helpful to understand the different abundance 
patterns between "strong" and "mild" Ba stars. 

There is a larger sample of Ba stars (318 from Lu cata- 
logue 1991) have been studied whose absolute magnitudes 
and kinematics by Gomez et al. (1997). Their results show 
that the samples of Ba stars are an inhomogeneous group, 
three of them among the four groups are disk population 
stars, and the fourth group contains halo stars. Mennessier 
et al. (1997) confirmed further this inhomogenity. If we can 
supplement the heavy-element abundances to compare the 
various characters of these different groups, it will be a 
stronger method to understand the natural properties of 
Ba stars. 

In addition, it is better to obtain the carbon and oxy- 
gen abundances together with s-process elements. The 
reason is that AGB stars are important nucleosynthesis 
sources of element carbon, which is taken out from stel- 
lar interiors together with s-process elements by the third 
dredge-up process (Liang et al. 2000; Liang et al. 2001). 
Element carbon will be then accreted by Ba stars, the 
companions of binary systems. 

Considering Ba stars from stellar population, it is in- 
teresting to combine Ba stars with CH stars. Some re- 
seachers suggest that metal-rich stars evolve into Ba stars 
and the metal-poor stars evolve to CH stars, which are 
metal-poor stars (Vanture 1992a, 1992b). At present, it is 
believable that nucleosynthesis of metal-poor AGB stars 
are more efficient than the metal-rich one (Busso et al. 
1999; Busso et al. 2001). Thus, it is valuable to study the 
heavy-element abundances of metal-poor CH stars. But 
CH stars are fainter, some are as low as 12 mag, which 
makes more difficult to get the required spectra. 

CH subgiants (Luck & Bond 1982, 1991; Smith et al. 
1993), Ba dwarfs (North et al. 1994) and extrinsic S stars 
(Jorissen et al. 1998) should be the closest groups to Ba 
giants. It will be interesting to combine and compare all 
of these group stars. Perhaps we will then be able to study 
Ba stars in an evolutionary sequence. 
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